Single-cell ionic current phenotyping explains stem cell-derived cardiomyocyte action potential morphology.

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) are a promising tool to study arrhythmia-related factors, but the variability of action potential (AP) recordings from these cells limits their use as an in vitro model. In this study, we use recently published brief (10 s), dynamic voltage-clamp (VC) data to provide mechanistic insights into the ionic currents contributing to AP heterogeneity; we call this approach rapid ionic current phenotyping (RICP). Features of this VC data were correlated to AP recordings from the same cells, and we used computational models to generate mechanistic insights into cellular heterogeneity. This analysis uncovered several interesting links between AP morphology and ionic current density: both L-type calcium and sodium currents contribute to upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the activation range of slow delayed rectifier K+ is the main determinant of AP duration. Our analysis also identified an outward current in several cells at 6 mV that is not reproduced by iPSC-CM mathematical models but contributes to determining AP duration. RICP can be used to explain how cell-to-cell variability in ionic currents gives rise to AP heterogeneity. Because of its brief duration (10 s) and ease of data interpretation, we recommend the use of RICP for single-cell patch-clamp experiments that include the acquisition of APs.NEW & NOTEWORTHY We present rapid ionic current phenotyping (RICP), a current quantification approach based on an optimized voltage-clamp protocol. The method captures a rich snapshot of the ionic current dynamics, providing quantitative information about multiple currents (e.g., ICa,L, IKr) in the same cell. The protocol helped to identify key ionic determinants of cellular action potential heterogeneity in iPSC-CMs. This included unexpected results, such as the critical role of IKr in establishing the maximum diastolic potential.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve.

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

Rapid ionic current phenotyping (RICP) identifies mechanistic underpinnings of iPSC-CM AP heterogeneity.

As a renewable, easily accessible, human-derived in vitro model, human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) are a promising tool for studying arrhythmia-related factors, including cardiotoxicity and congenital proarrhythmia risks. An oft-mentioned limitation of iPSC-CMs is the abundant cell-to-cell variability in recordings of their electrical activity. Here, we develop a new method, rapid ionic current phenotyping (RICP), that utilizes a short (10 s) voltage clamp protocol to quantify cell-to-cell heterogeneity in key ionic currents. We correlate these ionic current dynamics to action potential recordings from the same cells and produce mechanistic insights into cellular heterogeneity. We present evidence that the L-type calcium current is the main determinant of upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the excitable range of slow delayed rectifier K+ is the main determinant of action potential duration. We measure an unidentified outward current in several cells at 6 mV that is not recapitulated by iPSC-CM mathematical models but contributes to determining action potential duration. In this way, our study both quantifies cell-to-cell variability in membrane potential and ionic currents, and demonstrates how the ionic current variability gives rise to action potential heterogeneity. Based on these results, we argue that iPSC-CM heterogeneity should not be viewed simply as a problem to be solved but as a model system to understand the mechanistic underpinnings of cellular variability.